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Abstract 

We present a calculation for single-inclusive larg e-pr pion production in longitudinally 
polarized pp collisions in next-to-leading order QCD. We choose an approach where fully 
analytical expressions for the underlying partonic hard-scattering cross sections are ob¬ 
tained. We simultaneously rederive the corresponding corrections to unpolarized scatter¬ 
ing and confirm the results existing in the literature. Our results allow to calculate the 
double-spin asymmetry Af L for this process at next-to-leading order, which will soon be 
used at BNL-RHIC to measure the polarization of gluons in the nucleon. 


1 Introduction 


The measurement of the proton’s spin-dependent deep-inelastic structure function g\ by the 
EMC [1] more than a decade ago made once again the spin structure of the nucleon an exciting 
topic, which deservedly continues to spark large activity by both theorists and experimental¬ 
ists. The original result, that the total quark spin contribution to the nucleon spin is only of 
the order of about 20% has subsequently been confirmed by other experiments and is well- 
established today. For various reasons that we will not review here, gluons may very well play 
a more important role for the proton spin than quarks. Consequently, there is now a flurry of 
experimental activity aiming at measuring the polarization of gluons in the nucleon. In terms 
of a parton density, the required information is contained in [j|] 
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written in A + = 0 gauge, where x is the gluon’s light-cone momentum fraction of the proton 
momentum P + , and p F the factorization scale appearing in a hard process to which the gluon 
contributes. is the field strength tensor, and G^ v its dual. In more simple terms, A g(x, g F ) 
describes the difference in probabilities for finding a gluon with positive or negative hclicity in 
a proton with positive hclicity, at “resolution” scale fi F : 


Ag(x,g F ) = g p (x,g F ) -g±(x,p F ) , 
where superscripts (subscripts) denote the proton (gluon) hclicity. 


( 2 ) 


Deeply-inelastic scattering (DIS), Ip —> I'X , is a standard process for studying nucleon 
structure. However, it is not an ideal process for measuring the gluon content of the nucleon, 
due to the fact that the virtual photon in DIS couples directly only to quarks. Inclusive struc¬ 
ture functions therefore depend on the gluon density only through scale-evolution, and through 
higher orders in QCD perturbation theory. This explains why the existing polarized-DIS data 
have told us very little about A g |3|, f|. One may attempt to get access to the gluon density 
by selecting the photon-gluon-fusion process in DIS, which contributes to final states such as 
heavy flavor pairs, or high-transverse momentum (p F ) hadron pairs. Indeed, the COMPASS 
experiment p at CERN and Hermes f| at DESY follow this approach. Unfortunately, the 
rather low energy in these fixed-target experiments and the ensuing large systematic uncertain¬ 
ties in the theoretical predictions complicate these efforts significantly. Dedicated experiments 
at a possibly forthcoming future polarized ep collider, like the EIC [7|, would presumably make 
these channels more promising, however. 

The BNL Relativistic Heavy-Ion Collider RHIC || is able to run in a mode with polarized 
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protons. Very inelastic pp collisions will then open np unequaled possibilities to measure A g. 
RHIC has the advantage of operating at high energies (a/S = 200 and 500 GeV), where the 
theoretical description is under good control. In addition, it offers various different channels in 
which Ag can be studied, such as prompt photon production, jet production, creation of heavy 
flavor pairs, or inclusive-hadron production. In this way, RHIC is expected to provide the best 
source of information on Ag for a long time to come. 

The basic concept that underlies most of spin physics at RHIC is the factorization the¬ 
orem [£j, which states that large momentum-transfer reactions may be factorized into long¬ 
distance pieces that contain the desired information on the (spin) structure of the nucleon in 
terms of its parton densities such as Ag(x,p F ), and parts that are short-distance and describe 
the hard interactions of the partons. The two crucial points here are that on the one hand 
the long-distance contributions are universal, i.e., they are the same in any inelastic reaction 
under consideration, and that on the other hand the short-distance pieces depend only on the 
large scales related to the large momentum transfer in the overall reaction and, therefore, can 
be evaluated using QCD perturbation theory. The factorized structure forces one to introduce 
into the calculation a scale of the order of the hard scale in the reaction - but not specified 
further by the theory - that separates the short- and long-distance contributions. This scale is 
the factorization scale p F mentioned above. 

As an example, let us consider the spin-dependent cross section for the reaction pp —> nX, 
where the pion is at high transverse momentum, ensuring large momentum transfer. This is 
the reaction we study in the following. The spin-dependent differential cross section is defined 
as 

dAa = - [dc r ++ — da + ~~\ , (3) 

where again the superscripts denote the hclicities of the protons in the scattering. The statement 
of the factorization theorem is then: 

dAa = Y J dXa J dXb f dz c Af a (x a ,p F ) Af b (x b ,p F )D^(z c ,p' F ) 

CL^bjC 

x dA& c ab (x a P A , x b P B , Pn/Zc, HR, Hf, h’f) , (4) 


where the sum is over all contributing partonic channels a + b —> c + X, with dAa c ab the 
associated partonic cross section, defined in complete analogy with Eq. ([![), the hclicities now 
referring to partonic ones: 
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A few further comments are in order here. First, Eq. ([|) is actually a slight extension of 
the factorization theorem compared to what we stated above: the fact that we are observing a 
specific hadron in the reaction requires the introduction of additional long-distance functions, 
the parton-to-pion fragmentation functions D*. These functions have been determined with 
some accuracy by observing leading pions in e + e _ collisions and in DIS. Even though there is 
certainly room for improvement in our knowledge of the D™, we assume for this study that the 
fragmentation functions are sufficiently known. 

Secondly, we have displayed the full set of required scales in Eq. (||). Besides the factorization 
scale /ip for the initial-state partons, there is also a factorization scale for the absorption of 
long-distance effects into the fragmentation functions. In addition, we have a renormalization 
scale fiR associated with the running strong coupling constant a s . 

As mentioned above, the partonic cross sections may be evaluated in perturbation theory. 
Schematically, they can be expanded as 

dA Kb = dAd ab 0) + —' dAd ai 1] + • • • • ( 6 ) 

7r 

dAer^ 0 " 1 is the leading-order (LO) approximation to the partonic cross section and is, for our 
case of pion production, obtained from evaluating all basic 2^2 QCD scattering diagrams. 
It is therefore of order a^. The lowest order, however, can generally only serve to give a rough 
description of the reaction under study. It merely captures the main features, but does not 
usually provide a quantitative understanding. The first-order (“next-to-leading order” [NLO]) 
corrections are generally indispensable in order to arrive at a firmer theoretical prediction for 
hadronic cross sections. For instance, the dependence on the unphysical factorization and 
renormalization scales is expected to be much reduced when going to higher orders in the 
perturbative expansion. Only with knowledge of the NLO corrections can one reliably extract 
information on the parton distribution functions from the reaction. This is true, in particular, 
for spin-dependent cross sections, where both the polarized parton densities and the polarized 
partonic cross sections may have zeros in the kinematical regions of interest, near which the 
predictions at lowest order and the next order will show marked differences. 

There has been a lot of effort in recent years [10-15] to obtain the NLO corrections for the 
spin-dependent cross sections most relevant for the RHIC spin program. By now, essentially 
the only remaining uncalculated corrections are those for the partonic cross sections in Eq. (f|), 
i.e., inclusive pion production. These corrections will be presented in this paper. We emphasize 
that it is very appropriate to provide the NLO corrections at this time: it is planned for the 
coming RHIC run (early 2003) to attempt a first measurement of A g through exactly the spin 


3 


asymmetry 


An - 
n LL ~ 


dAa da ++ — da + 


(7) 

da da ++ + da + ~ 

for high-px pion production. The main underlying idea here is that A 1X LL is very sensitive to 
A g through the contributions from polarized quark-gluon and gluon-gluon scatterings. We 
note that the PHENIX collaboration has recently presented first, still preliminary, results for 
the unpolarized cross section for pp —> n°X at \fS = 200 GeV, which are well described by 


the NLO QCD calculation [16|, providing confidence that the theoretical framework based on 
perturbative-QCD hard scattering and summarized by Eq. ([4]) is adequate. 

Section 2 gives an outline of the calculation, summarizing the main ingredients. In Sec. 3 
we present some first numerical applications of our results. 


2 Calculation of the NLO corrections 

2.1 Outline of the strategy of the calculation 


The “parton-model” type picture employed in Eq. (|4|) implies that the partonic cross sections 
dAa^ b are single-inclusive cross sections for the reactions a + b —> c + X , i.e., summed over all 
final states (excluding c) possible at the order considered, and integrated over the entire phase 
space of X. Writing out Eq. (||) explicitly to NLO, we have 
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where Zq = 1 — V + VW, with hadron-level variables 

V = 1 +| , W = -0^ , S = ( Pa+Pb ) 2 , T = ( P A -P , U = (. P B -PV) 2 , ( 9 ) 

and corresponding partonic ones 
t — I! 

v = 1 + - , w = -— , S = (p a + Pb ) 2 , t=(p a - p c ) 2 , U = (p b - p c f . (10) 

s s + t 

Neglecting all masses, one has the relations 


s = x a x b S , t = —T , u = —U , x a = - , x h = —- -- . (11) 

vwz c z c { 1 - v) 
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The LO partonic cross sections dAa^°\v) are calculated from the 2 —> 2 QCD scattering 
processes, that is, X consists of only one parton, and its phase space is trivial and leads to the 
5(1 — w ) factor in Eq. (|j). We do not need to present the cross sections here, which have been 
known for a long time for both the unpolarized and the polarized cases O- There are actually 
only four generic reactions, qq' —> qq', qq —> qq , qq —> gg , and gg —> gg] all other processes 
follow from crossing if one works in terms of helicity amplitudes for each reaction, keeping all 
particles polarized. All tree-level 2 —» 2 helicity amplitudes are given in |L8|. The four generic 
processes give rise to the ten separate LO channels 


qq' —> qX 
qq' —> qX 
qq —> q'X 
qq —> qX 
qq —> qX 
qq -> gX 
qg -> qX 

qg gX 

gg -»■ gX 
99 qX , 


( 12 ) 


the “observed” final-state parton fragmenting into the hadron. At NLO, we have 0(a s ) cor¬ 
rections to the above reactions, and also the additional new processes 


qq' - 

- gX 

qq' - 

- gX 

qq - 

- gX 

qg - 

+ q'X 

qg - 

q'X 

qg - 

- qX 


(13) 


A singlc-inclusive-parton cross section is, of course, not a priori infrared-finite in QCD, but 
sensitive to long-distance dynamics through the presence of collinear singularities that arise 
when the momenta of partons in the initial or final states become parallel. Such a situation can 
appear for the first time at 0(al) (NLO), where 2^3 scattering diagrams contribute. From the 
factorization theorem discussed above it follows that long-distance sensitive contributions may 
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be factored into the bare parton distribution functions or fragmentation functions. The result of 
this procedure are finite partonic hard-scattering cross sections dAa^ b . At intermediate stages, 
however, the calculation will necessarily show singularities that represent the long-distance 
sensitivity. In addition, for those processes that are already present at LO, real 2 —> 3 and 
virtual one-loop 2 —> 2 diagrams contributing to the calculation will individually have infrared 
singularities that only cancel in their sum. Virtual diagrams will also produce ultraviolet poles 
that need to be removed by the renormalization of the strong coupling constant at a scale Hr. 
As a result, a regulator has to be introduced into the calculation that makes all the singularities 
manifest so that they can be canceled in the appropriate way. Our choice will be dimensional 
regularization, that is, the calculation will be performed in d — 4 — 2e space-time dimensions. 
Subtractions of singularities will generally be made in the MS scheme. 


Dimensional regularization becomes a somewhat subtle issue if polarizations of particles are 
taken into account. This is due to the fact that projections on helicities involve the Dirac matrix 
75 for quarks and the Levi-Civita tensor e pupa f or gluons. These two objects are genuinely four¬ 
dimensional and hence do not have a natural extension to 4 — 2e dimensions. In fact, some 
care has to be taken to avoid algebraic inconsistencies in the calculation when using 75 and 
e M " p<T . At the level of the algebra the treatment of 75 and e pupa of course only affects terms 
that are of 0(e). However, poles proportional to 1/e and 1/e 2 present in the calculation may 
combine with this to eventually result in non-vanishing contributions. A widely-used scheme 
for dealing with 75 and e pvpa in a fully consistent way is the one developed in |H|, the HVBM 
scheme. This is the scheme we have used for our calculation. It is mainly characterized by 
splitting the d -dimensional metric tensor into a four-dimensional and a (d — 4)-dimensional 
one. The Levi-Civita tensor is then defined by having components within the four-dimensional 
subspace only, and 75 anti-commutes with the other Dirac matrices in the four-dimensional 
subspace, but commutes with them in the (d — 4)-dimensional one. The HVBM scheme leads 
to a higher complexity of the algebra and of phase space integrals. However, one may make use 
of computer algebra programs such as Tracer 0 that allow to handle the split-up of space- 
time, and the treatment of (d — 4)-dimensional components in phase space integrals has become 
rather standard by now. We emphasize that for our present case the treatment of 75 and e pup(T 
has no bearing on the ultraviolet (renormalization) sector of the calculation, since we have no 
chiral vertices in the calculation. For instance, we may perform all renormalizations at the 
level of vertex and self-energy diagrams, without reference to the polarizations of the external 
particles. 


As remarked above, we need to sum over all possible final states in each channel ab —» cX, 
in compliance with the requirement of single-inclusiveness of the cross section. For instance, in 
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case of qg —> qX one needs, besides the virtual corrections to qg —> qg , three different 2 —> 3 
reactions: qg —> q(gg), qg —■> q(qq ), qg —■► q(q'q') (where brackets indicate the unobserved 
parton pair). Only all three processes combined will allow to arrive at a finite answer in the 
end. The summation over X is therefore always implicitly understood in the following. 


In addition, the two unobserved partons in the 2 —> 3 contributions need to be integrated 
over their entire phase space. The integration may be performed in basically two different 
ways. The first one relies on Monte-Carlo integration techniques. As was shown in |2T] , [ 22 || , the 
regions where the squared 2^3 matrix elements become singular can be straightforwardly 
identified and separated. These regions will yield all the poles in 1/e after integration, which 
eventually must cancel as described above. It then becomes possible to organize the calculation 
in such a way that the singularities are extracted and canceled by hand, while the remainder 
may be integrated numerically over phase space. This approach has the advantage of being 
very flexible; it may be used for any infrared-safe observable, with any experimental cut [^l] 


On the other hand, the numerical integration involved turns out to be rather delicate and 
time-consuming. In case of polarized collisions, the method was employed for the reactions 
pp —► jet A" [l3j and pp —► 7 A |L2| at NLO. 


The method we will employ is to perform the phase space integration of the 2 —> 3 contri¬ 
butions analytically. This has several advantages. In first place, the final answer is much more 
amenable to a numerical evaluation, giving much more stable results in a much shorter time. 
This may become important at a later stage, when experimental data will have been obtained 
and one is aiming to extract A g from them within a “global analysis” [^3j]. In addition, the 
“analytic method” has also been employed in the unpolarized case 0 . Since the calculation 
of the unpolarized and the polarized NLO terms largely proceeds along similar lines, we can 
compute both simultaneously. Our results for the unpolarized case may then be compared at 
an analytical level to those available in the numerical code of |24[ , which provides an extremely 
powerful check on the correctness of all our calculations. 


We will now separately address the virtual 2^2 and real-emission 2 —> 3 NLO contribu¬ 
tions. Then we will discuss their sum and the cancellation of singularities. 


2.2 Virtual contributions 

At O(af), virtual corrections only contribute through their interference with the Born diagrams, 
as sketched in Fig. [l]. We have calculated the virtual contributions with two different methods. 
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Figure 1: Interference of generic virtual (box, vertex, self-energy) contributions with Born diagrams. 


Firstly, we have performed a direct calculation. Here we could make use of known MS- 


renormalized one-loop vertex and self-energy structures as given in |25|, which may be readily 
inserted into the Born diagrams. One then additionally needs to calculate the box diagrams 
which are ultraviolet finite and hence not subject to the renormalization procedure. We have 
simultaneously computed the virtual corrections for the unpolarized case and found complete 


agreement with the results published in [26 


The second approach makes use of the fact that in Ref. [^7] the helicity amplitudes for 
all one-loop 2 —> 2 QCD scattering diagrams were presented. It is clear that these contain 
the information we need for our calculation. The only subtlety is that the helicity method 


employed in [g7| will not immediately yield the answer for the HVBM prescription we are 
looking for. However, as was also pointed out in [^, |27|, the translation between the results 
for the two schemes is fairly straightforward. In fact, by inspecting the singularity structure of 
the diagrams, one can derive a universal form for the virtual contribution V that schematically 
reads: 


V(s,t,u) = 


^ n n ) 


l 


m<n 


log — j B mn (s , t , u) + V(s, t, u ), 

(14) 

where n, m are summed over all external legs, the pi are the external parton momenta, and 
B denotes the Born cross section for the reaction under consideration. The B mn are the so- 
called “color-linked” Born cross sections, to be calculated according to rules given in |22j . The 
Ci and 7 j are coefficients depending only on the type of external leg, with C q = Cp = 4/3, 
C g = Ca = 3, 7 q = 3Cp/2, 7 g = /3 0 /2 = 11/2 — rif/ 3, rif being the number of flavors. Finally, 
V is the finite remainder. The only difference between the result for the virtual correction in 
the helicity amplitude method and the conventional HVBM scheme resides in the B and B mn 
terms. For the helicity method, these are four-dimensional quantities, whereas in conventional 
dimensional regularization they are calculated in d dimensions in the HVBM scheme. This 
property allows a direct determination of the full virtual correction in the HVBM scheme, since 


V has been calculated with helicity amplitude methods in [^]. This strategy for determining 
the virtual corrections was also adopted in Ref. 










































Figure 2: Representative 2 —>• 3 Feynman diagrams contributing to ab —► cA" to C>(a;^). 

We found complete agreement between the results obtained for our two approaches for 
obtaining the virtual corrections. 


2.3 2 —> 3 real contributions 


Figure [2| shows some representative 2^3 Feynman diagrams contributing to ab —> cA" to 
0(al). The squared spin-dependent matrix elements in d dimensions, using the HVBM pre¬ 
scription for 75 and the Levi-Civita tensor, are too lengthy to be reported here. Again, we 
have simultaneously calculated the squared matrix elements for the unpolarized case, and we 
recover the known [^ 6 | results in d dimensions. The polarized matrix elements can be checked 
in d — 4 dimensions against the expressions in [|H|], and again we find agreement. 
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In d — 4 — 2e dimensions, as a consequence of using the HVBM scheme with its distinction 
between four- and (d — 4)-dimensional subspaces, the squared matrix elements contain scalar 
products of vectors separately in these subspaces. For instance, while an outgoing unobserved 
parton with momentum k is massless, k 2 = 0, we may encounter its (d— 4)-dimensional invariant 
mass, denoted as k 2 , in the calculation, which is constructed from the (d — 4)-dimensional 
components of k. Such terms need to be carefully taken into account in the phase space 
integrations. 


The most economical way of organizing the phase space integration is to work in the rest 
frame of the two unobserved final-state partons whose momenta k 2 and k 3 can then be param¬ 
eterized as 


k 2 = (fc 0 , k 0 sin 6>i cos 0 2 ,k y , k 0 cos 6 U k) , 

k 3 = (ko, — ko sin 9\ cos # 2 , — k y , — ko cos$i, — k) , (15) 


and to define the momenta of the other three particles to lie in the x — z plane in the four¬ 
dimensional space. In this case the above k 2 is the only invariant arising from the (d — 4)- 


dimensional subspace. In Eq. (]T5|) k 0 = ^/s 23 /2 with s 23 = (k 2 + k 3 ) 2 = sv (1 — w), and k y 
denotes the unspecified y component of k 2 and k 3 which can be trivially integrated over since 
the matrix element does not depend on it. One then has the three-body phase space ||TT 


Ti — 


x 




(47T) 4 r(l - 2e) 

/*7T 

/ d0\ sin 1-2e 0\ I dd 2 sin 


dvv (l—v)~ 


dw (w(l — w)Y 


-2e 


dx 


—e) J 0 VT 


-X 


—(l+e) 


X 


(16) 


where x is k 2 normalized to its upper limit, x = Ak 2 /s 23 sin 2 9 1 sin 2 0 2 . 


The integrations we do analytically are over x (for those terms in the squared matrix ele¬ 
ments that have dependence on k 2 ) and the angles 6\ and 0 2 . v and w, defined in Eq. O. 
become integration variables in the convolution with the parton densities, according to Eqs. @ 
and ©• Extensive partial fractioning of the squared matrix elements always leads to the 
master integral for the angular integrations 


dB\ sin 1 2e Q\ I dd 2 sin Z£ 0 2 


-2 E, 


(1 — cos 9\)' J {1 — cos 9 1 cosx — sin 9 1 cos 9 2 sin x) 1 


T(1 — 2e) -_j5(l — e—j,l—e — l) nf.j, 2 X^ 

T(l-e) T(l-e) v 2 J ’ 


(17) 


where the last line is the result given in Ref. |28|. B(x,y) is the Euler Beta function and 
2 F 3 (a,b,c,z) denotes the hypergeometric function. 
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The final step in the evaluation of the 2 —> 3 contributions is to extract the poles arising 
when the invariant mass of the unobserved partons becomes small: S 23 —> 0. According to 
Eq. ( |T0| ) or the definition of .S 23 below Eq. ([L5D, this is the case for w —> 1. The fact that the 
LO contribution is proportional to 5(1 — w) indicates that the dependence on w is in the sense of 
a mathematical distribution. At NLO, the integrated matrix elements have terms proportional 
to 1/(1 — w) which, after inclusion of the factor (1 — w)~ £ from the phase space integral Q, 
can be expanded as 


^ ^-l-£ h,, \ 1 ( ln(l — w) 

(l-w) = --S(l- w) + Jr — f --e 


+ 0 (£ 2 ) . 


(18) 


\ 1 — w 

making the singularities at £ = 0 manifest. Here the “+”-distributions are defined in the usual 
way, 


f(w) [g(w)] + dw = / [f(w) -/(!)] g(w)dw . 


(19) 


2.4 Cancellation of singularities 


As mentioned earlier, genuine infrared singularities cancel in the sum of virtual and real con¬ 
tributions, among them all poles proportional to 1/e 2 which arise when soft and collinear 
singularities coalesce. The sum of virtual and real pieces is still singular for e —> 0 as a result 
of collinear divergencies. These remaining poles need to be factored into the bare parton dis¬ 
tribution functions and fragmentation functions, depending on whether their origin was in the 
initial or final state. 


Figure § sketches a typical collinear situation in a 2 —> 3 process. The contribution displayed 
will require a subtraction of the form oc 7 A P qq x A a qq ^ qq , where A P qq is the spin-dependent 


LO q —> q Altarclli-Parisi splitting function [£9| and Acr qq ^ qq represents the subsequent polar¬ 
ized LO scattering qq —> qq, evaluated in d dimensions. More precisely, the structure of this 
particular collinear subtraction is 


a 

7T 


— dx AH, 


m 




[xs, xt, u, e) 


dv 


5 (x (s + t) + u) 


( 20 ) 


where 

A H qq (z,fi F ) = + 7 e- In4 tt^ A P qq (z) (j^J + A /g<?(-) • ( 21 ) 

Here the Euler constant 7 e and In 47 T are the terms that are commonplace to subtract in order to 
work in the MS scheme. Af qq (z) is an additional finite piece in the subtraction that represents 
the freedom in choosing a factorization prescription and will be discussed below. We see in ([Li]) 
how the factorization scale fip emerges in the subtraction. In general, a process at NLO will 
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Figure 3: Representative collinear contribution to the subprocess qq —> qqg (see text). 


require several collinear subtractions, in both the initial and the final states. Depending on 
which types of partons are collinear, the other splitting functions A P qg , A P gq , A P gg , as well as 
other 2 —> 2 cross sections, will contribute. In the final-state collinear case, a singularity occurs 
when the observed parton and an unobserved one become collinear. The subtraction needed 
here can be easily written down in a form analogous to ( 2 T); it will involve the final-state 
factorization scale g! F . Note that, since we are not considering polarization in the final state, 
only spin-independent splitting functions appear in the final-state factorization subtraction. 

Taking the MS scheme literally, one would not have any additional finite pieces in the 
subtraction, beyond those involving and ln 47 r. That is, one would define (A )fij{z) = 0 
in the functions (A)//,_,■ involved in the various subtractions in the polarized and unpolarized 
cases. However, there is a well-known 11, 30, 31] subtlety arising in the q —> q splitting in 
the polarized case that is related to the use of the HVBM scheme. It is a property of the 
HVBM-scheme definition of 75 that it leads to helicity non-conservation at the qqg vertex in d 
dimensions. This can be seen from the non-vanishing difference of unpolarized and polarized 
d-dimensional LO quark-to-quark splitting functions: 


A Pt 2e {x) - PqZ 2e (x) = 4C F e (1 - x) 


m 


( 22 ) 


A disagreeable consequence of this is a non-zero first moment (rr-integral) of the flavor non¬ 
singlet NLO anomalous dimension for the evolution of spin-dependent quark densities, in con¬ 
flict with the conservation of flavor non-singlet axial currents p0| , [HI . |32|. It is therefore advis¬ 
able - albeit not mandatory in a purely mathematical sense - to slightly deviate from the MS 


scheme in the polarized case by choosing |30, 31 


Af qq (z) = -AC f (1-z) . (23) 

It is important to point out that the choice of the function A f qq (z) corresponds to the freedom 
in defining a factorization scheme. Of course, a physical quantity like the pion production 
cross section must not depend on the convention regarding which finite terms are subtracted 
from the partonic cross sections along with the collinear poles. Indeed, the parton distributions 
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functions are scheme dependent as well, so that at any given order in a s the scheme dependence 
cancels in the physical observable. The factorization scheme defined by the choice (|23|) has also 
been used in the available sets of spin-dependent NLO MS parton densities, so our definition is 
consistent with these densities. Since the HVBM “ 75 -effect” mentioned above is to be regarded 
as an artifact of the prescription and may be removed in a straightforward way by exploiting 
the conservation of non-singlet axial currents, results of polarized NLO calculations are usually 
regarded as being “genuinely” in the conventional MS scheme only with the choice (^3j). All 
other possible A f ig are, however, set to zero, as in the usual MS scheme. Needless to say that 
in the unpolarized case one has f qq = f qg = f gq = f gg = 0 in MS. 


2.5 Final results 

Once we have performed the factorization of collinear singularities, we arrive at the final result 
for the NLO partonic hard-scattering cross sections. We first of all note that, as mentioned 
earlier, we have calculated in parallel the NLO corrections for the unpolarized case. We have 
compared them term-by-term with the known analytical results in the code of |M| and found 
complete agreement. 

Our results for the spin-dependent NLO corrections may for each of the 16 subprocesses be 
cast into the following form: 


d ^ab l) { s ^ v ,w,n R ,n F ,n' F ) _ (a s (n R ) 


dvdw 


7T 
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—— f In — + A 2 5( 1 - w) In — 

(1 -w)+ J s s 
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(1 -w)+ V l ~ w J + 

+GTn(l — v) + H ln(l — w) + I ln(l — vw) + Jln(l — v + vw) 

In w 


+K 


1 — w 


L ln + M _l n ( 1 ~v + vw) 


1 — w 


1 — w 


(24) 


where all coefficients are functions of v and w, except those multiplying the distributions <5(1 — 
w), 1/(1 — w) +1 [ln(l — w)/(l — w)], which may be written as functions just of v. Terms with 
distributions are present only for those subprocesses that already contribute at the Born level, 
see Eq. (|I2|). 


We finally make a few observations about our results for the polarized case. Consider, 
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for example, the subprocess qq —> q'X in Eq. (p~2|). All Feynman diagrams contributing to 
this cross section at NLO, virtual as well as real, are annihilation diagrams, meaning that 
the initial quark and antiquark legs are part of the same fermion line. Independently of the 
number of gluons attaching to the fermion line, hclicity conservation in QCD demands that 
the annihilation can only occur if the quark and antiquark have opposite helicities. Keeping in 
mind the definition for the polarized cross section, we are led to the expectation that 




(25) 


should be fulfilled for this process. The only way in which this relation could be broken is if 
the regularization we adopt in the NLO calculation does not respect helicity conservation. As 
we discussed earlier, the HVBM prescription for y 5 indeed has this deficiency. However, as 
known from |3(], |3]J , the additional finite term A f qq (|23|) in the factorization subtraction (pT|) 
is precisely designed to cure this shortcoming of the HVBM prescription and to restore hclicity 
conservation. This is probably the most tangible reason why the choice Eq. (^3|) is required from 
a physical point of view. The implication of this is that our final results for qq —> q'X should 
indeed satisfy (|25|) , which we have verified. We can actually go one step further: the channels 
qq —> gX and qq —*■ qX have contributions from annihilation diagrams as well, but also ones 
from non-annihilation diagrams, for which the q and q scatter via f-channel gluon exchange. 
Helicity conservation makes no immediate statement about the non-annihilation diagrams. 
However, the channels qq' —>• gX and qq' —> qX are described by the non-annihilation diagrams 
alone. If we subtract the corresponding cross sections from the ones for qq —> gX and qq —> qX, 
respectively, we can use helicity conservation again for the remainder. Explicitly, we expect: 


dAa, 


9 ,( 1 ) 


99 


dAa. 


9,(1) 

99' 


da q - {1) 

u,u qq 


da q JP 


99 


(26) 


and similarly for an observed gluon. Again we have verified that our final results obey this 
relation, which we consider a very powerful check on the correctness of our results. 


3 Numerical results 

In this Section, we present a first numerical application of our analytical results. Instead 
of presenting a full-fledged phenomenological study of single-inclusive hadron production in 
polarized pp collisions, which we leave for a future study, we only report the main features of 
the NLO corrections and describe their impact on the cross sections and the spin asymmetry 
Predictions for A^ L are in immediate demand for an extraction of A g at RHIC in the 
very near future. 
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For our calculations we assume the same kinematic coverage as in the recent PHENIX mea¬ 
surement of the unpolarized cross section at a/S = 200 GeV [HJ, that is, we consider pion 
transverse momenta in the range 2 < px < 13 GeV and pseudorapidities \r )| < 0.38. We also 
take into account that the pion measurement is at present possible only over half the azimuthal 
angle. 


We will evaluate cross sections and spin asymmetries at both LO and NLO, in order to 
study the size and importance of the corrections we have calculated. We will always perform the 
NLO (LO) calculations using NLO (LO) parton distribution functions, fragmentation functions, 
and the two-loop (one-loop) expression for a s . To calculate the NLO (LO) unpolarized pion 
cross section needed for the denominator of the spin asymmetry A]° L in Eq. ( 0 ), we use the 
CTEQ5M (CTEQ5L) parton distribution functions. In all our calculations we use the 
pion fragmentation functions of Ref. |34|], which provides both a LO and a NLO set. For the 
polarized cross section, we will mainly use the (NLO/LO) “standard” sets of the spin-dependent 
GRSV H parton distributions (“GRSV-std”). Since we also want to investigate the sensitivity 
of A]° l to the polarized gluon density A g, we use another set of GRSV distributions, for which 
the gluon is assumed to be particularly large (“GRSV-max”). We note that the value of the 
strong coupling a s to be used in conjunction with the unpolarized parton distributions differs 
from that employed in the fits for the polarized sets and the fragmentation functions. Our 
convention will be to calculate the cross sections always with the strong coupling constant 
accompanying the parton distributions used. 

Figure [| shows our results for the unpolarized and polarized cross sections at NLO and LO, 
where we have chosen the scales /ir = /if = p'f — Pt■ The lower part of the figure displays the 
“It-factor” 

d(A)a NL ° 

d(A)a LO ‘ ( 7) 

One can see that in the unpolarized case the corrections are roughly constant and about 50% 
over the p t region considered. In the polarized case, we find generally smaller corrections which 
become of similar size as those for the unpolarized case only at the high-p^ end. The cross 
section for p t values smaller than about 2 GeV is outside the domain of perturbative calculations 
as indicated by rapidly increasing NLO corrections and, therefore, is not considered here. 


As we have mentioned in the Introduction, one reason why NLO corrections are generally 
important is that they should considerably reduce the dependence of the cross sections on the 
unphysical factorization and renormalization scales. In this sense, the A'-factor is actually 
a quantity of limited significance since it is likely to be rather scale dependent through the 
presence of the LO cross section in its denominator. The improvement in scale dependence when 
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Figure 4: Unpolarized and polarized 7 r° production cross sections in NLO (solid) and LO (dashed) 
at \/S = 200 GeV. The lower panel shows the ratios of the NLO and LO results in each case. 

going from LO to NLO is, therefore, a better measure of the impact of the NLO corrections, and, 
perhaps, provides also a rough estimate of the relevance of even higher order QCD corrections. 
Figure [] shows the scale dependence of the spin-dependent cross section at LO and NLO. In 
each case the shaded bands indicate the uncertainties from varying the unphysical scales in the 
range Pt/2 < hr = Pf = p'f < 2 px- The solid lines are for the choice where all scales are set 
to pt- One can see that the scale dependence indeed becomes much smaller at NLO. 

Finally, we consider the spin asymmetry which is the main quantity of interest here. Figure ^ 
shows A^ l , calculated at NLO (solid lines) for the “standard” set of GRSV parton distributions, 
and for the one with “maximal” gluon polarization We have again chosen all scales to be 
Pt- For comparison, we also show the LO result for the GRSV “standard” set (dashed line). 
As expected from the larger K factor for the unpolarized cross section shown in Fig. [|, the 
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Figure 5: Scale dependence of the polarized cross section for 7T° production at LO and NLO in the 
range p T / 2 < Pr = Pf = p'f — 2 Pt ■ We have rescaled the LO results by 0.1 to separate them 
better from the NLO ones. In each case the solid line corresponds to the choice where all scales are 
set to pr- 


asymmetry is somewhat smaller at NLO than at LO, showing that inclusion of NLO QCD 
corrections is rather important for the analysis of the data in terms of A g. 


We also conclude from the figure that there are excellent prospects for determining A g(x) 
from A]° l measurements at RHIC: the asymmetries found for the two different sets of polarized 
parton densities, which mainly differ in the gluon density, show marked differences, much larger 
than the expected statistical errors in the experiment, indicated in the figure. The latter may 
be estimated by the formula 


SA 


7T _ 

LL ~ 


l 


(28) 


where P is the polarization of one beam, C the integrated luminosity of the collisions, and Obm 
the unpolarized cross section integrated over the p^-bin for which the error is to be determined. 
We have used the very moderate values P = 0.4 and C = 7/pb, which are targets for the 
coming run. As mentioned above, we also take into account that at present with the Phenix 
experiment a pion measurement is possible only over half the azimuthal angle. 
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Figure 6: Spin asymmetry for 7T° production, using the “standard” set of GRSV @ and the one with 
“maximal” gluon polarization. The dashed line shows the asymmetry at LO for the GRSV “standard” 
set. The “error bars” indicate the expected statistical accuracy targeted for the upcoming run of 
RHIC (see text). 

4 Conclusions 

We have presented in this paper the complete NLO QCD corrections for the partonic hard- 
scattering cross sections relevant for the spin asymmetry A n LL for high -px pion production in 
hadron-hadron collisions. This asymmetry is a promising tool to determine the spin-dependent 
gluon density in the nucleon and will be measured in the coming run with polarized protons 
at RHIC. Our calculation is based on a largely analytical evaluation of the NLO partonic cross 
sections. 

We found that the NLO corrections to the polarized cross section are somewhat smaller 
for RHIC than those in the unpolarized case. The polarized cross section shows a significant 
reduction of scale dependence when going from LO to NLO. Upcoming RHIC data should be 
able to provide first information on A g even for rather moderate integrated luminosities. 
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Note added: 


While nearing completion of onr work, we learned that D. de Florian has performed [|35| the same 
calculation, using the “Monte-Carlo” method outlined in Sec. 2. This provides an extremely 
welcome opportunity for comparing the results. Early comparisons show very good agreement 
of the numerical results. 
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